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Spatial signal amplification in cell biology: a lattice-gas model for self- tuned phase ordering 
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Experiments show that the movement of eukaryotic cells is regulated by a process of phase separation of 
two competing enzymes on the cell membrane, that effectively amplifies shallow external gradients of chemical 
attractant. Notably, the cell is able to self-tune the final enzyme concentrations to an equilibrium state of phase 
coexistence, for a wide range of the average attractant concentration. We propose a simple lattice model in 
which, together with a short-range attraction between enzymes, a long-range repulsion naturally arises from 
physical considerations, that easily explains such observed behavior. 
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Specific moments of the life of a cell living in a multicel- 
lular community, such as migration, proliferation, organiza- 
tion in layers or complex tissues, imply spatial organization 
along some axis of direction. The original spatial symme- 
try of the cell must be broken to adapt to a highly structured 
anisotropic environment. For instance, migrating cells must 
orient towards sources of chemical attractants, mitotic cells 
must orient along the spindle-pole axis to bud daughter cells, 
epithelial cells must recognize the inner and outer part of tis- 
sues to define organ boundaries. From a physical point of 
view, these spatial organization phenomena may be seen as 
self-organized phase ordering processes, where the cell state, 
spontaneously, or because driven by an external field, decays 
into a state of coexistence of two or more chemical phases, 
spatially localized in different regions in order to define a front 
and a rear, a top and a bottom, an outer and an inner part. 
Local thermodynamic equilibrium requires precise tuning of 
chemical potentials to nongeneric values to allow the coex- 
istence of different phases. To implement this requirement 
in a robust way, capable of giving a stable response over a 
wide range of stimulation situations, biological systems must 
be endowed with self-organized tuning mechanisms leading 
to phase coexistence and polarization. 

Directional sensing in chemotacting eukaryotic cells pro- 
vides a beautiful illustration of these principles lll4ll . At the 
heart of directional sensing lies a chemical phase separation 
process taking place on the inner surface of the cell mem- 
brane O] . The main players of the process are the enzymes 
phosphatidylinositol 3 -kinase (PI3K), and phosphatase and 
tensin homolog (PTEN), which catalyze the switch of the 
phospholipid phosphatidylinositol between the bisphosphate 
(PIP2) and the trisphosphate (PIP3) states. The phospholipids 
are permanently bound to the inner face of the cell membrane, 
while the two enzymes diffuse in the cell volume and become 
active when they are adsorbed on the membrane. PI3K ad- 
sorption takes place through binding to receptors of the ex- 
ternal attractant. PTEN adsorption takes place through bind- 
ing to the PTEN product, PIP2, a process which introduces an 
amplification loop in the system dynamics iS [5|] . A second 
amplification loop provided by PI3K binding to PIP3 ||fl] has 
been recently observed. 



Although there are no relevant enzyme-enzyme or 
phospholipid-phospholipid interactions, the above described 
catalytic processes, together with phospholipid diffusion on 
the cell membrane, mediate an effective short-range attractive 
interaction among enzymes of the same type. This interaction 
drives the system towards phase separation in a PTEN-rich 
and a PI3K-rich phase ol |4J] . Two different regimes of mem- 
brane polarization may be distinguished. In the presence of 
an attractant gradient, anisotropy driven polarization is real- 
ized in a time of the order of a few minutes, and results in the 
formation of a PI3K-rich patch on the membrane side closer to 
the attractant source and of a PTEN-rich patch in the comple- 
mentary region ||5|] . The process works as an efficient gradient 
amplifier: a few percent gradient is sufficient to completely 
polarize the cell membrane. The orientation response is re- 
versible: by inverting the gradient direction the polarization 
orientation is also inverted. On the other hand, cells exposed 
to uniform distributions of attractant polarize in random di- 
rections over a longer timescale. The average concentration 
of attractant is of crucial importance, as shown by experimen- 
tally observed dose-response curves \m'. directional sensing 
does not take place neither at very low nor at very high attrac- 
tant levels, and there exists an optimal attractant concentration 
such that the cell response is maximal. 

On the basis of a simple analogy with the physics of bi- 
nary mixtures, one would expect that the coexistence between 
the PI3K-rich and the PTEN-rich phase would require a fine 
tuning of the chemical potential difference between the two 
species. Surprisingly, phase separation takes place instead 
for a wide range of absolute concentration of the attractant, 
and therefore of absolute values of the chemical potential for 
PI3K adsorption. To explain this mechanism we propose here 
a simple lattice-gas model in which, together with the effec- 
tive short-range attraction between enzymes, a long-range re- 
pulsion naturally arises from the finiteness of the enzymatic 
reservoir, that easily explains the observed behavior. 

Model - We represent the cell membrane by a square lat- 
tice of size L with N sites, using periodic boundary condi- 
tions. The sites / occupied by PI3K (PTEN) are described 
by a 5', = +1 (-1) spin Q. We denote by A^4 the total 
number of ± 1 enzymes in the cell, which is given by the sum 



of the number of cytosolic (free) enzymes and the number of 
membrane-bound enzymes: N^^ = A^^^^^ -^N^. The probabil- 
ity that a PI3K enzyme binds to site / is proportional to the 
number of cytosolic PBKs and to the density of binding sites 
(activated receptors with local concentration Ci and PIPs's). 
As a first approximation, the PIP3 concentration can be as- 
sumed to be linearly dependent from the density of PBKs. 
This gives, on site /: 



^(-1^+1)< 




A^free 



(1) 



where a+ , p+ , Cq" are functions of the chemical reaction rates, 
and di are the nearest neighbors of /. Similarly, the probability 
that a PTEN molecule binds to site / is proportional to the 
number of free PTENs, and to the concentration of PIP2: 



^(+1 -^ -l)oca 




free 



(2) 



We interpret A:^ =ln[^(-l ^+l)/^(+l -^ -1)] as an en- 
ergy difference (in units of RbT) between states St = -\-l and 
Si = —I, depending both on the local field Y^jedi ^j ^^^ ^^ the 
number of cytosoHc PBKs and PTENs. Since A^+ + A^" = A^, 
we can express N^^^ and A^^^^^ as functions of the magnetiza- 
tion m = {N^ —N~)/N. Linearizing A^ around Y^jedi^j = ^ 
and m = 0, we obtain: 



A^ = -2J Y, Sj - Ihi 

jedi 



2Xm 



where / 

ho 



1 / a+(3+ 

2 \ci^a+c+ 
r a-c-m+ \ 
V a+c+m- J ' 



and X 



= iln 

1 f 1 



+ ^) 



(3) 



-ho, with 



U-^ + -^), with 



2N^^/N — 1. If ^ < ^ we can neglect the dependence of 

J on the attractant concentration c/. 

Eq. ([3]) corresponds to the variation of the Hamiltonian 



^=-jY,SiSj- 

(iJ) 



l^htS^^-TSiSj. 



N 



(4) 



KJ 



The model dU contains a short-range ferromagnetic interac- 
tion representing the effective attractive interaction between 
enzymes, a long-range antiferromagnetic interaction which 
results from the finiteness of the cytosolic enzymatic reser- 
voir, and an external site-dependent field representing the ef- 
fect of the attractant. The latter depends on the concentra- 
tion Ci of activated receptors, which we take proportional 
to the concentration of external attractant, in the form Ci = 
c(l+£sin2Ssin2^). 

When hi is independent of /, the second and third term of 
Eq. dll) can be written (apart from a constant) as ^(^ — Jn)^, 
so that energy minimization leads the system to self-tune to 
the magnetization value h/X. Eq. (O shows that Si is sub- 
ject to the action of an effective external field /leff,/ = hi — Xm. 
The value /zgff,/ measures the degree of metastability of the 



PTEN phase, and tends to zero during the self-tuning evo- 
lution of the system. To realize the phase separation, / has 
to be greater than the critical value for the two-dimensional 
Ising spin model (/ :^ 0.44). Furthermore, in the absence 
of chemical attractant, the membrane has to be fully occu- 
pied by PTEN, implying ho > \. We set for definiteness 
J = ho = ^ = a+c^ = 1. 

Simulations - We study by Monte Carlo simulations the 
dynamics and the final state attained by the system, using a 
square lattice of size L = 2048. We first consider the case 
£ = 0, which corresponds to uniform stimulation. In the ab- 
sence of stimulation (c = 0, implying h = —I and m = — 1) the 
membrane is uniformly populated by PTEN molecules. Set- 
ting c > (which implies h > —1), spin up (PBK) domain 
nucleation is started in the spin down (PTEN) sea. The mag- 
netization m tends asymptotically to h, while the effective field 
/zeff tends to zero (Fig.[T]), realizing the condition for phase co- 
existence. 
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FIG. 1 : Self-tuning dynamics in the presence of a uniform activa- 
tion field h. The magnetization m grows to compensate the external 
activation field h. On the right, equilibrium states corresponding to 
different values of h. 



After a rapid nucleation phase, a domain coarsening dy- 
namics follows: large domains grow and smaller ones 
shrink fs^ (Fig. El). 
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FIG. 2: Coarsening dynamics leading to random membrane polar- 
ization in the presence of a uniform activation field. 

The final equilibrium state is characterized by the coexis- 
tence of the PBK and the PTEN phase, localized in two com- 
plementary clusters. The equilibrium position of the PBK 
cluster, which determines the direction of cell movement, 
is random. This behavior is consistent with experiments in 
which cells exposed to a uniform attractant distribution orient 
randomly (stochastic polarization) [T]. 

In the presence of a gradient in the chemical attractant 



(8 > 0) the PI3K cluster localizes around the maximum of the 
attractant density. To measure the polarization degree we de- 
fine the following order parameter: 



c> = - 



2 



-c)Si 



Ef|Q-c| 



(5) 



which is both a measure of the degree of order in the system 
and of the correlation of the center of the PI3K cluster with 
the maximum of the attractant density (Fig. [3]). 
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FIG. 3: Time evolution of the order parameter for different values 
of the activation field h, and for a fixed value 8 = 0.05 of the gra- 
dient. At the end of the polarization process the PI3K cluster (gray 
in the panels on the right) is centered around the point of maximum 
attractant stimulation (crosses). 



Dose-response curve - Simulations reproduce the qual- 
itative behavior of experimentally observed dose-response 
curves (tI |9|], showing no response for either very high or 
very low values of the attractant concentration, and optimal 
response for intermediate values (Fig.|4l). This effect can be 
explained as follows. For very low c the critical radius for 
patch nucleation is larger than the size of the cell, and no po- 
larization is possible. For very high c, such that /z > 1, the 
equilibrium magnetization is 1 , the whole system is uniformly 
populated by the PI3K phase, and again no polarization is pos- 
sible. Polarization is possible only for values which are inter- 
mediate between these two limit cases. 
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FIG. 4: a): Orientation degree of a population of cells as a func- 
tion of the attractant concentration, for a constant gradient (adapted 
from yji). b): Simulated equilibrium values of the order parameter a 
as a function of the attractant concentration c, for a constant gradient. 



Reversibility - Polarization induced by the gradient is re- 
versible. By changing the gradient direction after the sys- 
tem has reached equilibrium, the position of the PI3K clus- 
ter adjusts to the new direction in a finite time (not shown). 
This effect reproduces the observed reorientation of eukary- 
otic cells under varying attractant gradients observed in the 
experiments flO]. Interestingly, after changing the sign of the 
relative gradient we observed reorientation taking place by a 
collective movement of the PI3K cluster, and not by its evap- 
oration and successive recombination. 

Gradient amplification and polarization time - The tran- 
sient states are characterized by a coarsening dynamics with 
the appearance of scaling laws in the process of domain for- 
mation Jsi [ill |l2|] . Our simulations show that, for a condition 
of uniform distribution of attractant, in the initial coarsening 
stage the average cluster radius (r) grows approximately as 
t^^^ . In Fig. [5] the inverse length of the total cluster boundary 
is plotted against time llldl . 
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FIG. 5: Time evolution of the inverse length of the total cluster 
boundary for different values of the gradient 8. The dotted lines 
show the slope of the power-law behaviors characterizing the growth 
regimes dominated, respectively, by the uniform component of the 
attractant (^t^^^) and by the gradient (~ t). Arrows show the posi- 
tion of crossovers between the two scaling behaviors. 



We define the polarization time tp as the time for which the 
order parameter a reaches 90% of its equilibrium value. If 
the attractant is uniformly distributed the coarsening process 
stops when the average cluster radius becomes of the order of 
the cell size, r ^ L, implying that the spontaneous cell polar- 
ization time scales as tp ~ 1/L^. 

In the case of an attractant gradient we observe instead a 
double scaling behavior. For t < t^, where ^ is a crossover 
time depending on the amplitude of the gradient 8, cluster 
growth proceeds approximately as in the uniform case, while, 
for t > t^, the process of polarization becomes anisotropic, 
and the average cluster size grows approximately linearly in t 
(Fig. [5]). The presence of this double scaling law implies that 
the polarization time behaves as ^p ~ a + b/e + c/8^ (Fig.©. 
We can understand the double scaling law as follows. In the 
presence of an attractant gradient polarization takes place in 
two steps. In the initial (tuning) step the gradient of the attrac- 
tant is negligible with respect to the uniform component of the 
attractant and cluster growth is approximately unaffected by 
its presence. In the meanwhile, free enzymes shuttle from the 



cytosolic reservoir to the membrane, lowering the chemical 
potential for further cluster growth and effectively canceling 
out the effect of the uniform component of the attractant. This 
process continues until times of order t^, when only the effect 
of the gradient component is left. At this point, fast polariza- 
tion in the direction of the gradient takes place. 
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FIG. 6: Polarization time as a function of 8, and time evolution of the 
order parameter a for different values of 8 (inset). 

The anisotropic stage of cluster evolution leading to di- 
rected polarization occurs only if ^£ < tp. Otherwise, the pres- 
ence of a gradient of attractant becomes irrelevant and only the 
stage of isotropic patch growth actually occurs. The crossover 
time t^ increases with decreasing £ until it becomes of the or- 
der of tp, implying the existence of a lower threshold £th of 
detectable gradients. For 8 > £th anisotropy-induced polar- 
ization is much faster than spontaneous polarization. This 
explains the experimentally observed effect of gradient am- 
plification in chemotacting cells and the observation 12] of a 
lower threshold of detectable gradients, below which there is 
no directional sensing. Our results also confirm the theoretical 
predictions of 1I12I1 . 

Discussion - We have introduced a simple lattice-gas model 
for the process of eukaryotic directional sensing, which repro- 
duces important aspects of the observed phenomenology and 
sheds light on the underlying physical mechanism. The model 
maps signaling molecules and enzymes in spin variables, and 
the effective interaction between enzymes on the membrane 
into a ferromagnetic coupling. Enzymes shuttling from the 
cytosolic reservoir to the membrane is shown to provide a 
fundamental self-tuning mechanism which drives the system 
towards phase coexistence and polarization, by counteracting 
the effect of the external activation field. In the presence of 
an attractant gradient this mechanism cancels out the isotropic 
component of the attractant distribution in a first (tuning) stage 
of cluster growth, preparing the ground for fast directed polar- 
ization in the direction of the gradient in the next stage. The 
control provided by enzyme shuttling is encoded in the cou- 
pling of the effective magnetic field /zgff with the local order 
parameter m, thus realizing an effective long-range repulsion 
between enzymes and introducing in the model an element of 
self-organization. The existence of two distinct stages in clus- 



ter evolution when an attractant gradient is present is signaled 
in the model simulations by the emergence of a double power 
law for the time evolution of clusters of signaling molecules. 
This shows up in the dependence of directed polarization time 
from the gradient: for £ <C l,k scales as £~^. 

Previous models of eukaryotic polarization postulated the 
existence of a global inhibitor, which was needed to cancel 
out the average value of the attractant leaving only the gradi- 
ent component (see lll3n for a review of previously proposed 
models of eukaryotic chemotaxis). The weak point of this ap- 
proach is the necessity to fine-tune the activity of the global 
inhibitor in order to attain perfect cancellation of the average 
attractant value. In our scheme instead, the exchange of chem- 
ical factors between the cell membrane and a finite cytoso- 
lic reservoir realizes a self-tuning mechanism which naturally 
leads to an equilibrium state characterized by the coexistence 
of two distinct phases, similarly to what happens in the case 
of first-order phase transitions in a closed liquid-gas system 
or in the precipitation of a supersaturated solution. The sens- 
ing mechanism encoded in the model is particularly robust, 
allowing the cell to respond over a wide range of attractant 
concentrations. 

Our model explains the experimentally observed behavior 
of chemotacting cells and reproduces several effects, such 
as gradient amplification, typical dose-response curves, re- 
versibility of orientation. More generally, it shows that im- 
portant biological functions may be described at a physical 
level as self-organized phase transition processes. 

We gladly acknowledge useful discussions with S. Di Talia, 
I. Kolokolov, V. Lebedev and G. Serini. 



* Electronic address: teresa.ferraro@na.infn.it| 
[1] A. Ridley etal. Science 302, 5651 (2003). 
[2] D.A. Lauffenburger and A.E Horwitz, Cell 84, 359 (1996). 
[3] A. Gamba, et al Proc. Natl. Acad. Sci U.S.A. 102, 16927 

(2005). 
[4] A. de Candia et al, Sci. STKE 378, pU (2007). 
[5] M. lijima and R Devreotes, Cell 109, 599 (2002). 
[6] M. Dance etal, J. Biol. Chem. 281, 23285 (2006). 
[7] S.H. Zigmond, J. Cell. Biol. 75, 606 (1977). 
[8] E.M. Lifshitz, and L.R Pitaevskii, Physical Kinetics (Pergamon 

Press, 1981). 
[9] L. Song et al, Eur. J. Cell Biol. 85, 981 (2006). 
[10] C. Janetopoulos etal, Proc. Natl. Acad. Sci. U.S.A. 101, 16606 

(2004). 
[11] A. Bray, Adv. Phys. 45, 357 (1994). 
[12] A. Gamba et al, Phys. Rev. Lett. (2007). 
[13] R Devreotes and C. Janetopoulos, J. Biol. Chem. 278, 20445 

(2003). 
[14] The biological facts presented below are taken, where not oth- 
erwise specified, from the reviews ^ and |2l. 
[15] One can imagine performing a coarse-graining of the system on 
an appropriate length scale and associating a +1 sign to PI3K- 
rich sites and a —1 sign to PTEN-rich sites. 
[16] If the system is composed of circular domains, the inverse 
length of the total cluster boundary scales as the mean radius 
of the clusters. 



